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Abstract 

. Lax pairs with operator valued coefficients, which are explicitly connected by means 

of an additional condition, are considered. This condition is proved to be covariant 
with respect to the Darboux transformation of a general form. Nonlinear equations 
arising from the compatibility condition of the Lax pairs in the matrix case include, in 
particular, Nahm equations, Volterra, Bogoyavlenskii and Toda lattices. The examples 
of another one-, two- and multi-field lattice equations are also presented. 
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I Introduction 



In the present paper we consider a large class of the integrable systems of nonlinear equations 
taking values on an associative ring of noncommutative operators. They are defined as 
the compatibility condition of Lax pairs characterized by the property that the equation 
for the time evolution of the wave function is explicitly determined (in a local way) by 
the coefficients of the spectral problem. There exist different types of Lax pairs, whose 
coefficients are connected explicitly. Lax pairs with coefficients defined on an associative 
algebra of scalar pseudo-differential operators were introduced in [p], || (see also || and 
references therein). These Lax pairs play an important role in the Sato theory [|J and in 
constructing the modifications of KP hierarchy [pi |(J. The case of the shift operators and 
associated lattice equations were discussed in 0. The approach, which is applied below to 
connect the coefficients of the Lax pairs, differs from known ones. The hierarchy of Darboux 
covariant nonlinear "multi-field" equations we describe in this paper contains, as the simplest 
case, the "one-field" equations of von Neumann type 

ip=[H{p),p], ip=[HJ(p)}. (1) 

These equations and their solutions were investigated in the context of the density matrices 
and Hamiltonians in || |9|, [10], where, for instance, the formulas of the Darboux transforma- 
tion were constructed. In the matrix case, the multi-field equations we derive here admit 
the reductions leading to known and new integrable nonlinear lattice systems. 

The technique exploited in this paper combines and develops the approaches of works 
TT, [T2[]. In Sec. II we show that the relations between the coefficients of the equations forming 



the Lax pair, i.e. necessary conditions for the Lax pair compatibility, are identically satisfied 
if an additional condition on the coefficients is imposed. This condition connects explicitly 
the coefficients of the equations of the Lax pair and allows us to write in a closed form 
the nonlinear equations, which follow the compatibility condition. A theorem establishing 
the Darboux covariance of the Lax pair with the additional condition and, consequently, of 
the corresponding nonlinear equations is proved in Sec. III. This theorem gives an effective 
tool of producing the infinite hierarchies of solutions, including the multi-soliton solutions, 
for nonlinear equations and their reductions. The equations of von Neumann type and 
associated lattice equations are considered in Sec. IV. The next two sections are devoted to 
two- and multi-field generalizations of the von Neumann type equations. The examples 
presented include known lattice equations as well as some new ones. Particular cases of the 
Darboux transformation satisfying the conditions of the theorem of Sec. Ill are discussed in 
the Appendix. 



II Lax pair and mult i— field equations 

Let us begin with the overdetermined system of linear equations (Lax pair) 

-iip = t/jA(X) , v 

z\ip = i>H{\) ■ { ] 

Here A and z\ are complex numbers, ip takes values in a given linear space £, A(X) and H(X) 
are linear operators £ i— > £ belonging to an associative ring, the dot denotes a derivative 
(i.e. an operator satisfying the Leibnitz rule). The compatibility condition of Eqs. (Q) is 

iH(X) = [A(X),H(X)}. (3) 
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If we assume that the operators entering the Lax pair are rational functions of A with operator 
valued coefficients of the form 

L M 1 

A{\) = £A fc £ fc + £-C7,, (4) 

k=0 k=l 
N 

H(X) = "£X k H k . (5) 

fc=0 

then Eq. ([3D becomes equivalent to the following system of algebraic and differential relations 
between operators B k , C k and H k : 

N 

£ [B m - k ,H k ] = (N <m<L + N), (6) 

fc=max{0, m—L} 
min{N,m+M} 

£ [C k . m ,H k ] = (-M<m<0), (7) 

k=0 

m min{A r ,m+M} 

iH m = E [B m - h ,H h ]+ £ [C fc _ m ,# fc ] (0<m<iV). (8) 

fc=max{0, m—L} fc=m+l 

The connection between operators B k , C k and ilfc, which is implied by Eqs. (HEO), is 
implicit. It is possible to express B k and C k explicitly in terms of the operators H k in order 
to satisfy Eqs. (j^,^) identically. Indeed, let us put 



1 / jL-k 

1 ' " -' NtjI-1\ -1^ 



1 / d M ~ k 

Ck = w^-[d^ 9{H{£),£) 



(9) 

<r=0 

(10) 



e=0 



where f(X, A) and A) are well defined functions of an operator X and parameter A. We 
assume hereafter that for any operator X(X), which is analytic in the neighborhood of the 
point A = oo, the function f(X(X), A) is also analytic and the condition 

[f(X(X),X),X(X)] = 0. (11) 

is valid in this neighborhood as well. In the case of the function g(X, A), analogous properties 
have to take place in the neighborhood of the point A = 0. In particular we have 

[g(X(X),X),X(X)] = 0. (12) 

Eqs. (||,0) are fulfilled for any B k and C k defined by Eqs. fl9|,|i~0|) as a consequence of the 
following identities 

jN+L-m 

lj ^ 1 —\f^H{<r x \ r 1 ), S N H{,-% =Q = (N<m<L + N), 

jM+m 

[g(H(e),e),H(e)] n = (-M < m < 0). 



Using Eqs. (|^,0) we can rewrite Eqs. as 



N m 

iH m = E [B m -k, H k ] + E[Cfe-m, Hk], (13) 

k=m+l k=0 



where the coefficients B k and for k < are calculated accordingly to Eqs. (p|JT0|) . The 
operator A(X) given by (§J^,0) is conveniently represented in the following equivalent way 

A(X) = [F(H(X), X)U + [G(H(X), A)] . (14) 

Here 

F(H(X),X) = X L f(H(X)/X N ,X), (15) 

G(H(X),X) = \- M g(H(X),X)i (16) 

symbols [...]oo and [...]o denote the parts of the power expansions in A that contain nonneg- 
ative and negative powers respectively. These equations establish a connection between the 
coefficients of Lax pair (|). 

In what follows we restrict our consideration to functions f(X,X) and g(X, A), which 
possess an additional property, namely they are covariant with respect to a similarity trans- 
formation T: 

f{T- x XT, A) = T-7(X, A)T, g^XT, A) = T^g^X, X)T. (17) 



Conditions (|IT ,112,1?) are not very restrictive and are satisfied, for example, by polynomials 
in X and sums of negative powers of polynomials in X with scalar coefficients. If X is 
selfadjoint operator, then these conditions are valid for all f(X, A) and g(X, A) determined 
via the spectral theorem. 




Ill Darboux covariance 

Let us consider the transformation 

m=^PD(X), (18) 

where ...[1] denotes the image under the transformation, and D(X) is an invertible linear 
operator depending on A. We say that the Lax pair (|2|) is Darboux covariant with respect to 
transformation {ip, A(X), H(X)} — > {ip[l], A[1](X), H [1](A)} if the following equations hold 

= mmw [ } 

and the structure of A[1](A) and if[l](A) is the same as the structure of A(X) and H(X). 
The notion of "structure" means that the shapes of the coefficients of Lax pairs (0) and 
(|il]|) are the same. The most important point is that the locations of singularities of A(X) 
and A[1](A), H(X) and if [1](A) and their types should coincide. The transformations of the 
form ( |il|[20| , |21"D that satisfy these conditions are called Darboux transformations . These 
transformations allow one to generate the hierarchies of solutions of nonlinear equations 
admitting the compatibility condition representation and of associated Lax pairs. In finite 
dimensional (matrix) cases D(X) is termed the Darboux matrix [14 . 

Substituting (|T8|) into flT9| ) we obtain expressions for the coefficients of the transformed 
Lax pair 

A[1](X) = -iD(X)- l D(X) + D(X)- l A(X)D(X). (20) 

H[1](X) = D(X)- 1 H(X)D(X). (21) 

If -D(A) and .D(A) -1 are regular on the plane of parameter A at singular points of the coeffi- 
cients of Lax pair then the sufficient condition for D(X) to define the Darboux transformation 
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comes from the requirement that the right hand sides in Eqs. ( p0| , |21|) have no the singulari- 
ties at the points, which are the singular point of D(X) and _D _1 (A). We refer the Reader to 
the Appendix, where the examples of such transformation are presented. 

The following theorem gives sufficient conditions of covariance of Eq. (|i~4|), which explic- 
itly connects the coefficients of the Lax pair (fj), with respect to the Darboux transformation 
(0,0,0). 

Theorem. If D(X) and D(\) 1 are rational functions in A that have poles at finite points 
/ii, . . . , n s and . . . , Us respectively, Hk 7^ (k = 1, . . . , S) and [D(A) _1 -D(A)] 00 = 0, then 
Eq. $1^) is Darboux covariant. 

Proof. Substituting (Mj into (M) yields 



A[1](A) = -iD{\y l b + D(\)- l [F(H(\), A)] ooJ D(A) + D{X)- l [G{H{X), X)] D(X). (22) 
It is enough to show that 

A[1](A) = [F(H[1](X), A)]oo + [G(H[1}(\), A)] . (23) 



The right hand side of ( P2|) is a rational function of A with poles at most at fii,...,ng 
and 0, oo. Therefore the following decomposition holds 

41] (A) = X>[l](A)k + [A[1](X)U + [A[l](X)]o (24) 
fe=i 

where [v4[l](A)] Alfc is the principal part of the Laurent expansion of A[1](A) at the point 
A = /ifc. Since the first term in formula (|22|) can have poles at most at /ii,...,/i5 and 
vanishes as A — > oo, we obtain 

[A[l](A)]oo = [D(A)- 1 [F(i/(A), A)] 0OJ D(A)] 0O + [D(A)- 1 [G(iJ(A), A)] £>(A)U 

[A[1](A)] = [D(X)- 1 [F(H(X)^)] CO D(X)} Q + [D(X)- 1 [G(H(X),X)} D(X)} . 
Using equalities 

[,D(A)- 1 [F(if(A),A)] 00J D(A)] 00 = [D(A)" 1 F(ff(A), X)D(X)} 00 , 
[ J D(A)- 1 [F(i/(A), A)] oo D(A)] = [[ J D(A)- 1 F(i/(A), A) J D(A)] oo ] = 0, 
[ J D(A)- 1 [G(iJ(A), A)] 0J D(A)]oo = [[ J D(A)- 1 G(if(A), X)D(X)] ] oo = 0, 
[D(A)- 1 [G(iJ(A), A)] D(A)]o = ^(A^G^A), A)D(A)] . 

we have 

[A[l](A)]oo = [,D(A)- 1 F(if(A), A) J D(A)] oc , (25) 

[A[1](A)] = [D(A)- 1 G(i/(A), A)D(A)] . (26) 
It follows from the definition of the Darboux transformation that for k = 1, . . . , S 

[A[l]{\)] n = 0. (27) 
Combining (|23|-f27p and taking into account Q2TJ) , (|TJ|-[I7|), we get (^). ■ 
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It is well known that the conditions ( |27| ) can be solved resulting in an explicit expression 
for D(\) in terms of solutions of the Lax pair (§) and a dual pair, which belong to the kernel 
of this operator or its inverse one. For this reason, we prefer to use the name of the Darboux 
transformation technique instead of the dressing method. In the Appendix we give examples 
of the Darboux transformations, which satisfy the conditions of our Theorem and can be used 
to construct the hierarchies of solutions of nonlinear equations (|j) under constraints (|9| JT0|) . 
The applications of the Darboux transformation technique to certain nonlinear equations of 
von Neumann type, including an equation in infinite dimensional case, which are of interest 
in connection with quantum mechanics and statistical physics, can be found in |8|, p|, [10 . 



If operator T in Eqs. ( |TTD is independent of A, then the expressions in right hand sides 
of these formulas correspond to so-called gauge transformation of wave function 

ip — > tp = ip T. 

The case of the Lax pairs with g(X, A) = is gauge equivalent to the case f(X, A) = if T 
solves equation 

if = B T. 

Some of the lattice equations presented in the next sections are the compatibility condition 
of gauge equivalent Lax pairs. 



IV Equations of von Neumann type (one— field equa- 
tions) 

In this section we consider the case N — 1, i.e., 

H(X) = XHt + H . 
The compatibility condition of the Lax pair gives us the equation 

iH Q = [B Q ,H ] + [C 1 ,H 1 ], (28) 



where 



c=o 



Ci 



1 



d 



M-l 



e=0 



(29) 



(30) 



(M-l)! \de M - 

We refer to this case as "one-field" because Hi is the constant operator. The right-hand side 
of Eq. (|28| ) combines both types of the nonlinearities as in the equations of von Neumann 
type (|l|). The nonlinear equations corresponding to simplest choices of functions f(X, A) and 
g(X, A) are presented below. We also obtain the equations that follow them if the matrix 
coefficients Hi and Ho are defined in a special manner. 



IV.l f(X, A) = iX l (/ G N), g(X, A) = 0, L = 1 
The compatibility condition leads to the equation 

Ho=j2[Hr- 1 H H l 1 - m ,H }. 

m=l 



(31) 
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This equation with I = 2 is multi-dimensional Euler's top equation [|15| , [T6| , I7[ . Darboux 
covariance of Eq. fl3~if ) and associated Lax pair was proved in || . 
Let matrices Hi and ifo have the form 



(32) 



Then Eq. ([H]) yields 



2-i 



Pk — Pk^2 (Pk+m — Pk~ 



(33) 



m=l 



These equations are known as the Bogoyavlenskii lattice ||18|| . In the case I = 2 Eqs. 
coincide with the Volterra system [Tj]| 



Pk — Pk{Pk+l — Pk-l), 



(34) 



which describes stimulated scattering of plasma oscillations by ions [ED 



IV.2 f(X, A) = iX~ n (n G N), g(X, A) = 0, L = 1 
In this case Eq. (|2~8D is written as given 



tt \^\H- m H Tjm-n-1 tt 

H — —2_^[ n i -no-ni 



o • 



(35) 



m=l 



The Bogoyavlenskii lattice (|33| ) with Z = n + 1 follows this equation if matrices Hi and H 
are chosen in the form 

Hl t kj = $k,j-l, H 0: kj = PkSk,j-n-l- (36) 

IV.3 /(X, A) = 0, g{X, A) = iX l (I e N), M = 1 

From Eq. (^8|) we have 

H =[H l ,Hi]. (37) 
Assuming that matrices Hi and ifo are represented in the following manner 

Hi,kj = fik,j-i+ii Ho,kj — Pk8k,j+i, (38) 



we obtain the well known lattice [18 



2-1 2-1 

Pk = H Pk-m — Y\. Pk+m- 
m=0 m=0 



(39) 



These equations with I = 2 are obviously reduced to the Volterra system ([34]). If we put 

p k = exp(u fc ) , 

then Eqs. fl3~9|) read as 

2-1 2-1 

u k = exp( Uk-m) - exp( u k+m ). (40) 

771=1 717=1 

As it was noted at the end of Sec. Ill this case is gauge equivalent to the case IV. 1. 
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IV.4 f(X, A) = 0, g(X, A) = iX l (I e N), M = 2 
Eq. ([28D yields 

Supposing 



m=0 



(41) 
(42) 
(43) 

m=0 i=0 i=0 i=0 i=0 

(It is assumed hereafter that 11^0 • ••i = 1 if w < 0.) In the case I = 3 these equations are 
equivalent to the Volterra system 



we have 



TO— 1 



l-m-2 



l—l m— 1 l—m—2 
Pk = J! ( IT Pfc-2i IT Pk+2i-l+2 — IT Pk-2i+l-2 H Pfe+2) 



IV.5 /(X, A) = 0, #(X, A) = (/ G N), M = 1 

Eq. fl28|) takes the form 

If matrices H\ and Hq are defined as follows 

Hi t kj — fik,j+i+i, Ho t kj = PkOk,j+u 

then the compatibility condition implies 

i i 

Pk — H Pk+m ~ IT Pk-m- 
m=l m=l 

Introducing two sets of new dependent variables 

p k = exp(-w fc ), 
Vk = p k \ 

we obtain equivalent representations of Eqs. 



(44) 
(45) 

(46) 



i i 
ii k = exp( Uk- m ) - exp( ^ u k+m ) 

m=0 m=0 



(47) 



and 



= V k( IT ~ IT V k+m)- (48) 

m,=l m=l 

Eqs. (|47D with Z = 2 were studied in [21]. Lax pair for Eqs. (48) was found in |18|. These 
equations in the case I = 1 obey a symmetry v k — » — ffc and look like a natural generalization 
of the Volterra system ([34]). 
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IV.6 f{X, A) = 0, g(X, A) = iX~ l (I E N), M = 2 
In this case Eq. (j28|) is written in the next manner 

/ 



Ho=~J2 [H^^K' 1 - 1 , H,]. (49) 

m=l 



If matrices ifi and H are defined as follows 

Hl,kj = <5fcJ+«+2, ^O.fcj = Pk$k,j+2, (50) 

then we come to equations 

l m l—m m l—m 

h = s ( n /^i;-^ n Pfc-2i- 2 - n pkUi n Pfc-w*) • (51) 

m=l i=l j=0 i=l i=0 

IV.7 f(X,\) = Q, g(X,\)=g(X), M=l 
Here we have 

i H = [g(H ),H 1 }. (52) 

By construction, g(Ho) commutes with Hq. The Lax representation and Darboux covariance 
properties of Eq. ( |52] ) with arbitrary well-defined function g(X) were established in [PH[ . 



The cases <7pT) = iX 3 and g(X) = iX 1 were considered in |22| in the framework of the 



symmetry approach to the classification problem of integrable equations on free associative 
rings. The Lax representation for the equations in these cases seems to be new. 

IV.8 f(X J X) = X\ g(X 1 X) = 0, L = 2 

The compatibility condition (fffi ) becomes 

iH = [h(H Q ),H ) = [H 1 ,F(H )) (53) 

(compare with Eqs. (M), where 

h{H ) = HqH^ + HqHxHqHx + H HfH + HiH 2 H 1 + HiHqH x Hq + H^H^, 

F(H ) = H Q Hi + HqHiH + HqHxHI + H 1 Hq. 
Let us note that, contrary to the previous example, [F(Hq),Hq] 7^ 0. We refer to the maps 



Hq 1 — > F(Hq) of such a kind as nonabelian functions , or nonabelian nonlinearities [|1^]. This 
example is a particular case of the equations ((4~T|) . 



V Two— Field Equations 

A few examples of systems appearing if N = 2 are considered in this section. 
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V.l f(X,X)=iX,g(X,X) = 0,L = l 

The compatibility condition (||) leads to equations 

Hq = [Hi,H ], 
Hi = [H 2 ,H ], 
H 2 = 0. 

It is checked immediately that functions 

Fx = (H 2 -aH )/(2i), 
F 2 = (H 2 + aH )/2, 
F 3 = H 1 /(2i), 

where a is a parameter (cr / 0), satisfy equations 

Fi = [F^Fzl+iiF^Fx], (54) 
F 2 = [F 3 ,F 1 ]+i[F 3 ,F 2 ], (55) 
F 3 = a[F u F 2 }. (56) 

In terms of new dependent variables 

fk = gF k g- 1 (A: = 1,2,3), 

where g solves equation 

g = -i g F 3 , 

Eqs. ([54 H561) are rewritten as 



fi = L/2,,/3], 

f2 = [/3 > AL 

fs = v[h,h\- 

If we impose condition f£ = — then a has to be real. This system with cr = 1 is known 
as Nahm equations J23], It will be shown in the next subsection that this case is also 
connected with Toda lattice equation. 

V.2 /(X, A) = 1X 1 (I e N), g(X, A) = 0, L = 1 

From Eqs. (Q) we have 

H = ^[H^HxH^^Hoi (57) 

m=l 

H, = [H l 2 , H ] + £ \E™- x H x H l 2 m , H x ], (58) 

m=l 

H 2 = 0. 



If we put 



H2,kj — Sk,j-i, H\ y kj — hkSkj+i-i, H 0t kj — Pk8kj+2i-i, (59) 
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then Eqs. (|57| , [58| ) read as 



1-1 



hk — Pk+i — Pk + h k ^2 {hk+m — hk-m), 



i-i 



m=l 



(60) 



Pk — Pk/j (hk+ m — hk-l-m)- 



m=0 



Let 



hk — &k- 



If coefficients pk are chosen as given 

(C is arbitrary constant), then system ( pDj ) is equivalent to the following equations 

a k = C(e^m=o( a k+m+l - Ok-m) _ e EL=o( (T fc+"i _ Vk-l-m)\ _|_ ^ ^ ((Xfc+m ~ ^k-m 



m=l 



Assuming I = 1 and C = 1 we come to the Toda lattice equation [^5], ^6], ^7 



(61) 



(62) 



Eqs. ( plf ) can be viewed as a generalization of the Toda lattice on the systems of particles 
interacting with finite number of nearest neighborhoods. For n x n matrices, these equations 
admit additional reductions 

<y m +i+k = -cr m +i-k if n = 2m + 1, 

a m+k = -a m+ x-k if n = 2m 

or 

a m+k = -<7 m -k if n = 2m 
(k = 0, ...,m). In the case I = 1, these reductions lead to generalized periodic Toda lattices, 



whose connection with the root systems of semisimple Lie algebras was established in p8 



V.3 f(X, A) = iX~ l (I 6 N), g(X, A) = 0, L = 1 

The compatibility condition in this case yields 

i 

tt \ 1 f TJ~ m U ITrn—l—l it 1 

11 — — 2^ l l li 2 n l n 2 i U 0\, 
m=l 



zj [rr-i tt 1 \ f TT—m tt Tjm—l—1 tt 1 

m=l 

H 2 = 0. 

If H2, Hi and Hq are defined in the following manner 

H^kj — $k,j+i, Hi,kj = hk5kj+i+i, Ho t kj = Pk$k,j+2i+i, 



(63) 
(64) 



(65) 



11 



then Eqs. fl63| , |64| ) are written as 

i 

hk = Pk+l — Pk — hk (hk+m — hk-m), 
m=l 

I 

Pk = Pk J! (hk-l-m ~ hk+rn)- 



(66) 



m=l 



In the case I = 1 these equations are so-called Belov-Chaltikian lattice [25|. The bilinear 
approach was applied to Belov-Chaltikian lattice in |30 . 



Expressing dependent variables in the terms of new ones 

hk = Ofc, 
p k = C r eSL=i( cr fc-i-m - °k+m) 
(C is a constant) we reduce Eqs. (|66"D to the Toda-type lattice equations 

a k = C (e^™=l( ak - m ~ a k+l+rn) _ e T,' m =i( a k~l-m ~ Cfc+m)^) _ ^ ^ (&k+m ~ &k- 



(67) 



m=l 



V.4 f(X, A) = 0, g(X, A) = iX l (I e N), M = 1 

In this case Eqs. (||) are rewritten as given 

ifn = [H l ,Hi], 



Hi 
H 2 



[<H 2 ] 
0. 



(68) 
(69) 



If matrices H 2 , H 1 and H have the form 

H 2 ,kj — Ok t j-2l+li Hl,kj — hkSkj-l+l, 

then Eqs. (MM) yield 



i-i i-i 

hk = J\ Pk-m — J\ Pk+l+m, 

m=0 m=0 

l-l l-l 



H 



0,kj 



pkS, 



fcj+i) 



(70) 



(71) 



Pk — hk~l Y\_ Pk-m — hk Y\_ Pk+r. 



m=0 



m=0 



These equations with / = 1 are equivalent to the Toda lattice (|62[) . 
V.5 f{X, A) = 0, g(X, A) = iX l (I e N), M = 2 

From Eqs. (|8|) we have 

z-i 

#o = [H , H 2 ) + [H™H 1 H l ~ m ~ 1 , Hi], 



m=0 



/-l 



#1 

H 2 



Y.lKHiH 1 ^- 1 ,^] 

m=0 

0. 



(72) 
(73) 
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Taking matrices H 2 , Hi and H as given 

H2,kj — Ok,j-2l+2> H\,kj — hkSk,j-l+2i ^O.fcj = Pk$k,j+2i 

we put compatibility condition into the form 



(— 1 m— 1 l—m—2 m—l l—m—2 

hk = X! ( hk-2m Y\. Pk-2i Y\. Pk+2i-l+2 ~ hk-2m+l-2 J\ Pk-2i+l-2 

m=0 i=0 i=0 i=0 /-H 

I— 1 m—l £— m— 2 

Pfc = ^fc-i J! hk-2m Y\. Pk-2i Y\. Pk+2i-l+2 ~ 
m=0 i=0 i=0 

2-1 m-1 Z-m-2 J-l J-l 

— hk X! ^/«-2m-2 17 Pk-2i-2 Pk+2i-l + JJ Pfc-2i — 17 Pfc+2i- 

m=0 t=0 i=0 i=0 i=0 



In the case I = 2 and hk = these equations coincide with the Volterra system (|3~4 ) 



V.6 /(X, A) = 0, flf(X, A) = (/ G N), M = 1 

In this case Eqs. yield 

i?i = [# ~ , H 2 ] , 
# 2 = 0. 

Let the matrices H 2 , Hi and H be represented in the following manner 

^2,kj — 3k,j+2l+li Hl,kj — hkOk,j+l+l, Ho,kj — Pk"k,j+l- 

The compatibility condition leads to the system 

i i 

hk = Yl Pk+m ~ II Pk-l-mi 



m=l 



m=l 



p k - hk+i n p k i +m - h k n Pk-m- 



m=l 



m=l 



V.7 f(X, A) = 0, g(X, A) = iX~ l (I 6 N), M = 2 

Eqs. (|) give 

i 

tt frr-i tj 1 \ f Tj—m tt Tjm—l—1 tt 

H o — [-"o )-"2j— 2^1-"o -"i-"o >■": 

m=l 

#! = -j2[H G - m H 1 H^- l - 1 ,H 2 ], 



m=l 



Supposing 



Ho = 0. 



— $k,j+2l+2, Hl,kj — hkSkJ+l+2, Ho } kj — Pk8k,j+2, 
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we have 

l m l—m m l—m 



- Y ( hk+2m~2l-2 II Pk+2i-2l-2 II Pk-2i-l-2 ~ h k +2m JJ Pk+2i II Pk-2i+l / 
m=l i=l i=0 i=l i=0 

I m l—m 

Pk = hY h k+2m -l-2 II Pk+2i-l-2 II Pk-2i-2 ~ ( 83 ) 



m=l i=l i=0 

I m l—m I 



- hk+i y h k+2m n p k+2i n p k \ i+ i + n p k \ 2i - n Pk-u- 

m=l i=l i=0 i=l i=l 

V.8 f{X, A) = iX 2 , g(X, A) = 0, L = 2 

The compatibility condition (^) in this case is written as 

H = [Hi, Hq] + [H 2 , H 2 ], 
Hi = [H 2 , H Hi + H\Hq], 
H 2 = 0. 

Taking matrices H 2 , Hi and H as follows 

H 2 ,kj — &k,j-ii Hi k j = h k 5 k j, H ok j = p k S k j + i, (84) 

we obtain 



h k — p k +\{h k+ i + h k ) — p k {h k + /ifc-i), 
Pk = Pk(pk+i - Pfc-i +h 2 k - h 2 k _ x ). 



J5) 



This system is the first member of the hierarchy of higher Toda lattices [[U], |32|]. In the case 
h k = this system is nothing but the Volterra system (|34T). 

VI Multi-Field Equations 

In this section we present nonlinear equations that follow Eq. (^|) with positive integer N for 
special choices of functions f(X, A) and g(X, A). 

VI.l f{X, A) = %X l (I E N), g(X, A) = 0, L = 1 
The compatibility condition is 

i 

Hi = [H l N , H^] + Y Wn~ 1 Hn- i N l ^ m , Hi] (86) 

m=l 

(i = 0,...,N). If we put 

Hi ,kj — Pk $kj+(N-i)l-l (87) 

(Pk = 1), then Eqs. ( |S6D give 



.(i) _ _ (t) ^ / (JV-l) _ (JV-l) v 

PA: — Pk+l Pk < Pk Z^\Pk+m Pk+{i-N+l)l-m) ' 

m=0 
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VI.2 f(X, A) = iX~ l (I 6 N), g(X, A) = 0, L = 1 
In this case Eqs. are written in the following way 

Hi = [H„ l , H^] - £ [HJ^Hn-iN^- 1 - 1 , Hi] (89) 
(i = 0, N). Assuming 

H itk j = pi'Skj-w-ty-i, Pk = ( 90 ) 

we obtain 

.(i) _ (*-i) (i-i) , (i) >p _ JJV-i)n / q n 

FA; — Pk-l Pk i FA; vFA:+(Af-i-l)«+m Pk—m )■ V ,L ) 

m=l 

VI.3 /(X, A) = 0, A) = (/ 6 N), M = 1 

The compatibility condition yield 

H t = [H l ,H i+1 ] (92) 
(z = 0, iV). Let matrices ifj have the form 

Hi ,kj — Pk ^kJ-U+1, (93) 
where p k = 1. In this case Eqs. (|92]) lead to the following lattice equations 

Pk =Pk-i n Pk—m Pk n Pk+ 

(94) 

m=0 m=0 

The case I = 1 was studied in [[33|, [34 



VI.4 f(X, A) = 0, g(X, A) = iX" 1 (/ e N), M = 1 
From Eqs. (|8|) we have 

Hi = [H-\Hi +1 ] (95) 
(z = 0, N). If matrices i7j are defined as follows 

Hi t kj = Pk $k,j+il+li (96) 

then Eqs. fl95|) give 

Pfc = Pk+l 1 1 Pfc+m — Pk II P*-«-n»> (97) 

m=l m=l 

where we use the notation 

(0) 

Pk = Pk- 

VII Conclusion 

In future we will continue the study of the von Neumann type equations and their gener- 
alizations presented in the previous sections. The integrals and the multi-soliton solutions 
will be considered. An investigation of the hierarchies of symmetries and compatible flows 
for these equations can lead to new hierarchies of integrable equations |22|, |35|, |36], |37f . The 
results will be of special interest in the case of integrable lattice equations. It should be also 
mentioned that discretizations of the lattice equations attract large attention in recent years 
(see, e.g., [BBf and references therein). 
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Appendix. Darboux Transformations 

Here we discuss briefly some particular cases of the Darboux transformation, which satisfy 
the conditions of Theorem in Sec. II. Let the operator D(X) in (|TB| ) be represented as given 

D(X) = 1 + ^-P, (98) 
li — A 

where 



X is a solution of the Lax pair (0) with parameter v: 

f ~ix = xA(u) 
\ ZvX = xH(u) ' 

(p is a solution of the dual Lax pair with parameter \i: 

f itp = A{(j,)<p 
1 = H(p)<p ' 

(x, ip) is a scalar product. It is obvious that P 2 = P and 

-iP = PA(u)P ± - P ± A(/j,)P 

(Pj_ = 1 — P). If coefficients of the transformed Lax pair flT5| ) are defined by 

L M i 

A(A)[1] = E^W + Eti^W. (99) 

fc=0 fc=l A 

2V 

#(A)[1] = EA^[1]> 

fc=0 



where 



fl fc [l] = B fc + (//-!/) E (v m ~ k - 1 P±B m P-v m - k - 1 PB m P ± ) 

m=k+l 
M 

C k [l] = C k -(fiL-v) E {jt^PlCnP ~ V k - m - X PC m P^) , 

V 

= F & + (//-*/) E (v m - k - 1 P±H m P-v m - k - 1 PH m P 



m=k+l 



then Eqs. ([T9p are identically fulfilled. This statement can be proved by direct computation. 
Since 

D(\r = i+^p, 

v — A 
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operators D(\) and -D(A) -1 have poles in points \x and v. It is seen from Eq. ( p9[ ) that Eq. 
(^) are valid. 

The formulas written above form the so-called binary Darboux transformation. The 
corresponding Darboux transformation for an n-dimensional matrix case is produced from 
them if we assume in (^8|) that 

p = ^(x^r 1 *, 

where \ and <p are respectively m x n and n x m matrix solutions of direct and dual Lax 
pairs. Some examples of Darboux transformations in infinite dimensional case, which are 
suitable for integrable lattice equations, were given in |13| . Very recently a new construction 
of the Darboux transformation in terms of Clifford numbers was described in [^] . 
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